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ABSTRACT 

In this paper we present an improved theoretical scenario concerning near infrared and 
visual magnitudes of RR Lyrae variables, as based on up-to-date pulsating models. 
New relations connecting V and K absolute magnitudes with periods, mass, lumi- 
nosity and metal content are discussed separately for fundamental and first overtone 
pulsators. We also show that the V — K colors are predicted to supply tight constraints 
on the pulsator intrinsic luminosity. On this basis, we revisit the case of the prototype 
variable RR Lyr, showing that the parallax inferred by this new pulsational approach 
appears in close agreement with HST absolute parallax. Moreover, available K and 
V measurements for field and cluster RR Lyrae variables with known reddening and 
metal content are used to derive a relation connecting the K absolute magnitude to 
period and metallicity (AfR--[Fe/H]-logP) as well as a new calibration of the My- 
[Fe/H] relation. The comparison between theoretical prescriptions and observations 
suggests that RR Lyrae stars in the field and in Galactic Globular Clusters (GGCs) 
should have quite similar evolutionary histories. The comparison between theory and 
observations also discloses a general agreement that supports the reliability of current 
pulsational scenario. On the contrary current empirical absolute magnitudes based on 
the Baade-Wesselink (BW) method suggest relations with a zero-point that is fainter 
than predicted by pulsation models, together with a milder metallicity dependence. 
However, preliminary results based on a new calibration of the BW method provided 
by Cacciari et al. (2000) for RR Cet and SW And appear in a much better agreement 
with the pulsational predictions. 
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1 INTRODUCTION 

Cosmic distances are fundamental cornerstones in modern 
astronomy, and reliable distance indicators can supply tight 
constraints on several cosmological parameters (see e.g. Goo- 
bar et al. 2000; Ferrarese et al. 2000). Accordingly, the num- 
ber of theoretical and empirical investigations looking for 
new standard candles and/or aimed at improving the classi- 
cal ones is countless in the recent literature (see e.g. Salaris 
& Cassisi 1997; Bono et al. 1999; Feast 1999; Walker 1999; 
Caputo et al. 2000; Carretta et al. 2000; Castellani et al. 
2000; Gieren, Fouque, & Storm 2000; Layden 2002). How- 
ever, the ongoing problems affecting distance determinations 
are soundly revealed by current uncertainties on the true 
distance modulus fj,o of the nearest galaxy, the Large Mag- 
ellanic Cloud (LMC), which is often taken as a reference for 



extragalactic distance determinations. In fact, current esti- 
mates range from ho,lmc ~ 18.37 ± 0.23 mag (RR Lyrae, 
Luri et al. 1998) to ~ 18.7 mag (classical Cepheids, Feast & 
Catchpole 1997; Bono et al. 1999; Caputo et al. 2002). 

In this paper we discuss RR Lyrae stars, the radial pul- 
sators connected with low-mass He-burning stars in the Hor- 
izontal Branch (HB) evolutionary phase. The use of these 
variables as distance indicators appears still hampered by 
significant uncertainties, and indeed the relation between 
absolute visual magnitude and iron content, My-[Fe/H], 
calibrated according to different empirical or theoretical 
constraints, supplies absolute magnitudes that for a given 
metallicity cover « 0.3 mag. Oddly enough, one also finds 
that the empirical approach based on the BW method, gives 
different results when different bands (optical vs near in- 
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frared) or different approaches (surface brightness vs canon- 
ical BW) are used. 

In this context, the endeavor to link RR Lyrae pulsation 
properties, i.e. periods and/or amplitudes of the light curve, 
to the absolute magnitudes of these stars (see Sandage 1993; 
Caputo 1997; Caputo et al. 2000; Caputo, DegPInnocenti 
& Marconi 2001) appear of great relevance, because peri- 
ods (and amplitudes) are firm and safe observational pa- 
rameters, independent of distance and reddening. Empirical 
investigations (Longmore, Fernley, & Jameson 1986; Long- 
more et al. 1990) have already suggested that in the K 
band the RR Lyrae stars in GGCs follow a rather well- 
defined Period-Luminosity (PLk) relation whose slope is, 
within observational errors, almost constant when moving 
from metal-poor to metal-rich clusters. This scenario has 
been recently supported by the theoretical investigation by 
Bono et al. (2001, Paper I). In that paper, on the basis of 
nonlinear, convective RR Lyrae models, we found that the 
if-magnitude of fundamental (F) and fundamentalized first- 
overtone (FO) pulsators can be fitted by a common Period- 
Luminosity-Metallicity (PLZk) relation which appears to 
be in good agreement with observations. 

In this paper we present improved predictions concern- 
ing both F and FO pulsators. Section 2 presents new theoret- 
ical relations concerning V and K absolute magnitudes. On 
this basis, in the same section we discuss improved PLZk 
relations, showing that V — K colors might allow us to con- 
strain the intrinsic bolometric luminosity of the pulsators. 
Section 3 deals with the prototype variable RR Lyr, used as 
fundamental calibrator of the RR Lyrae distance scale, while 
observational data for field and cluster RR Lyrae stars are 
discussed in Section 4. In Section 5 we present a new relation 
connecting the K absolute magnitude to period and metal- 
licity (Mx-[Fe/H]-logP) as well as a refined My-[Fe/H] rela- 
tion. A few final remarks concerning both theoretical models 
and observations are outlined in Section 6 



2 NEW PLZk AND PLCvk RELATIONS FOR F 
AND FO PULSATORS 

To properly implement the pulsational framework presented 
in Paper I, we computed additional pulsating models adopt- 
ing a finer grid of effective temperatures. We ended up with 
132 F and 91 FO nonlinear models that cover a wide metal- 
licity range 0.0001 < Z < 0.02. For each metal abundance 
the helium content was fixed according to a primordial he- 
lium abundance of F p =0.24 and a helium-to-metal enrich- 
ment ratio AY/AZ « 2.5. 

At fixed chemical composition, we explored three dif- 
ferent luminosity levels that cover current uncertainties on 
the ZAHB luminosity as well as post-ZAHB evolutionary ef- 
fects (see column 3 in Table 1). As for as the pulsator mass 
is concerned, we adopted for each given metallicity a "ref- 
erence mass" M r , roughly following the correlation between 
metal content and mass of ZAHB models populating the RR 
Lyrae instability strip, i.e. 3.77 < logT e < 3.86. These M r 
values are listed in Table 1 (column 4), together with the 
evolutionary mass (with an average uncertainty of ~ 2%) of 
ZAHB models at logT e =3.85 and 3.80, as inferred according 
to evolutionary computations available in the literature. In 
order to investigate the dependence of the predicted rela- 



Table 1. Input parameters for the grid of RR Lyrae models. Mass 
and luminosity are in solar units. 



z 


Y 


log L 


M r 


M(3.85) 


M(3.80) 










±2% 


±2% 


0.0001 


0.24 


1.61, 1.72, 1.81 


0.75 


0.796 


0.852 


0.0004 


0.24 


1.61, 1.72, 1.81 


0.70 


0.699 


0.721 


0.001 


0.24 


1.51, 1.61, 1.72 


0.65 


0.648 


0.666 


0.006 


0.255 


1.55, 1.65, 1.75 


0.58 


0.585 


0.589 


0.01 


0.255 


1.51, 1.57, 1.65 


0.58 


0.575 


0.578 


0.02 


0.28 


1.41, 1.51, 1.61 


0.53 


0.545 


0.546 



tions on the pulsator mass, additional models were finally 
constructed at selected metal contents (Z = 0.0001 and 
Z = 0.001) but with slightly different reference mass val- 
ues. 

The entire set of fundamental and first overtone non- 
linear models were integrated in time till they approach 
the limit cycle stability. These calculations provided peri- 
ods, pulsation amplitudes, and bolometric light curves. The 
light curves were transformed into the observational plane 
by adopting bolometric corrections and color-temperature 
relations provided by Castelli, Gratton & Kurucz (1997a,b). 
Tables 2 and 3 list the predicted periods (in days) and mean 
(magnitude-weighted) V and K magnitudes for the new F 
and FO pulsation models, respectively. These new data im- 
plement similar predictions already presented in Tables 2 
and 3 of Paper I. 

Data listed in Table 2 and Table 3 disclose that, for 
any given mass, luminosity and metal content, the visual 
magnitude My of the pulsators is quite independent of 
the effective temperature, whereas Mk shows a strong de- 
pendence on this parameter (see also Paper I). This ef- 
fect is the consequence of the fact that in the range of 
RR Lyrae effective temperatures the visual bolometric cor- 
rection is roughly constant, whereas in the if-band it is 
strongly dependent on the effective temperature. It follows 
that the if-band magnitudes become significantly brighter 
when moving toward cooler effective temperatures. The pul- 
sation equation P = f(L,M,T e ) connects, for each given 
chemical composition, periods to luminosity, mass and ef- 
fective temperature. Therefore, the dependence of the ab- 
solute magnitude on T e (through the BC) can be replaced 
with a dependence on period. Accordingly, when moving 
from the blue to the red edge of the instability strip the pe- 
riod increases at roughly constant absolute visual magnitude 
(5Mv/8\ogP ~ —0.14), whereas Mk shows a much steeper 
correlation (8M K /8\ogP ~ -2.1.) 

On the other hand, when varying the luminosity for 
any given value of mass and effective temperature, one has 
again a relation between absolute magnitude and period. 
The key feature of the K-band is that the two derivatives 
(5MK/8\ogP) M ,L) and {8Mk /8\ogP) M ,T c ) are quite simi- 
lar. In other words, one finds a sort of degeneracy for which 
the near-infrared period-magnitude [PLk) relation is only 
marginally dependent on the bolometric luminosity of pul- 
sators (see Fig. 1 in Paper I). 

As for the dependence of the PLk relation on the pul- 
sator mass and metal content, current calculations confirm 
the results presented in Paper I. However, at variance with 
the approach adopted in that paper, where F and funda- 
mentalised FO periods (i.e., log Pf = log Pfo + 0.127) were 
taken simultaneously into account, we found that the intrin- 
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Table 2. Periods and absolute mean (magnitude- weighted) mag- 
nitudes for the new fundamental models. 



log L/L c . 



(K) 



P 

(d) 



M K 
(mag) 



M v 
(mag) 



0.0001 


1 


.61 


6800 


0.3876 


-0.172 


0.848 


0.0001 


1 


.61 


6700 


0.4077 


-0.226 


0.845 


0.0001 


1 


.61 


6600 


0.4291 


-0.282 


0.844 


0.0001 


1 


.61 


6500 


0.4518 


-0.328 


0.834 


0.0001 


1 


.61 


6400 


0.4764 


-0.378 


0.833 


0.0001 


1 


.61 


6200 


5299 


-0.470 


0.843 


0.0001 


1 


.61 


6100 


0.5594 


-0.507 


0.848 


0.0001 


1 


.72 


6800 


0.4790 


-0.448 


0.597 


0.0001 


1 


.72 


6600 


0.5292 


-0.554 


0.572 


0.0001 


1 


.72 


6400 


0.5882 


-0.658 


0.562 


0.0001 


1 


.72 


6200 


0.6561 


-0.752 


0.566 


0.0001 


1 


.81 


6600 


0.6304 


-0.778 


0.357 


0.0001 


1 


.81 


6400 


0.7005 


-0.880 


0.338 


0.0001 


1 


.81 


6300 


0.7399 


-0.930 


0.335 


0.0004 


1 


.81 


6600 


0.6632 


-0.778 


0.330 


0.0004 


1 


.81 


6500 


0.7014 


-0.830 


0.320 


0.0004 


1 


.81 


6400 


0.7400 


-0.881 


0.315 


0.006 


1 


.75 


6900 


0.5988 


-0.494 


0.428 


0.006 


1 


.75 


6800 


0.6290 


-0.546 


0.426 


0.006 


1 


.75 


6700 


0.6614 


-0.598 


0.416 


0.006 


1 


.75 


6600 


0.6972 


-0.650 


0.108 


0.006 


1 


.75 


6500 


0.7349 


-0.702 


0.402 


0.006 


1 


.75 


6400 


0.7755 


-0.753 


0.399 


0.006 


1 


.75 


6300 


0.8184 


-0.803 


0.400 


0.006 


1 


.75 


6200 


0.8626 


-0.850 


0.403 


0.006 


1 


.75 


6000 


0.9632 


-0.962 


0.404 


0.006 


1 


.75 


5900 


1.0184 


-0.988 


0.423 


0.006 


1 


.75 


5800 


1.0794 


-1.029 


0.430 


0.01 


1 


.65 


6800 


0.5162 


-0.300 


0.697 


0.01 


1 


.65 


6600 


0.5708 


-0.403 


0.682 


0.01 


1 


.65 


6500 


0.6019 


-0.455 


0.672 


0.01 


1 


.65 


6400 


0.6322 


-0.509 


0.668 


0.01 


1 


.65 


6300 


0.6669 


-0.558 


0.665 


0.01 


1 


.65 


6200 


0.7052 


-0.607 


0.664 


0.01 


1 


.65 


6100 


0.7441 


-0.658 


0.668 


0.01 


1 


.65 


5900 


0.8345 


-0.762 


0.681 


0.01 


1 


.65 


5700 


0.9331 


-0.839 


0.695 



sic accuracy can be improved by independently treating the 
two pulsation modes. 

By accounting also for the mass dependence, a linear 
regression through the models gives: 



Table 3. Periods and absolute mean (magnitude- weighted) mag- 
nitudes 

for the new first overtone models. 



z 


logL/L Q 




P 


M K 


My 






(K) 


(d) 


(mag) 


(mag) 


0.0001 


1.61 


7100 


0.2546 


-0.031 


0.817 


0.0001 


1.61 


7000 


0.2657 


-0.081 


0.817 


0.0001 


1.61 


6900 


0.2784 


-0.131 


0.812 


0.0001 


1.61 


6800 


0.2915 


-0.180 


0.808 


0.0001 


1.61 


6600 


0.2935 


-0.276 


0.806 


0.006 


1.75 


7000 


0.4132 


-0.453 


0.325 


0.006 


1.75 


6900 


0.4349 


-0.502 


0.350 


0.006 


1.75 


6800 


0.4563 


-0.552 


0.363 


0.006 


1.75 


6700 


0.4806 


-0.603 


0.370 


0.006 


1.75 


6600 


0.5070 


-0.656 


0.374 


0.006 


1.75 


6500 


0.5343 


-0.709 


0.378 


0.006 


1.75 


6400 


0.5615 


-0.759 


0.383 


0.006 


1.75 


6300 


0.5904 


-0.806 


0.389 


0.006 


1.75 


6200 


0.6222 


-0.849 


0.394 


0.01 


1.51 


6800 


0.2871 


0.045 


0.965 


0.01 


1.57 


6900 


0.3070 


-0.055 


0.800 


0.01 


1.57 


6800 


0.3228 


-0.104 


0.807 


0.01 


1.65 


7000 


0.3394 


-0.204 


0.574 


0.01 


1.65 


6900 


0.3562 


-0.253 


0.594 


0.01 


1.65 


6800 


0.3745 


-0.302 


0.603 


0.01 


1.65 


6700 


0.3931 


-0.355 


0.611 


0.01 


1.65 


6600 


0.4135 


-0.402 


0.615 



Ml 



0.511-2.102 log P+0.095 log Z-0.734 log L- 



1.753 log 
(1) 



M F K ° 



0.029-2.265 log P+0.087 log Z-0.635 logL-1.633 

(2) 

where mass and luminosity are in solar units and the rms 
scatter is a=0.016 mag. As already discussed in Paper I, 
the predicted slope at constant luminosity and metallicity 
appears in good agreement with empirical data for GGC 
variables, as given, e.g., by the RR Lyrae stars in M3 col- 
lected by Longmore et al. (1990). 

The present enlarged set of pulsating models confirms 
that over the explored range of metallicity and luminos- 
ity the Mk — log P correlation is only mildly dependent on 
the luminosity, a variation of 0.1 dex implying a change of 



8Mk ~ 0.07 mag, as well as on the metallicity, a varia- 
tion of 0.3 dex implying a change of SM K ~ 0.03 mag. The 
predicted magnitudes are thus expected to be only mildly 
dependent on uncertainties introduced by these theoretical 
or observational parameters. On the other hand, the non 
negligible dependence on the pulsator mass should only in- 
troduce a marginal uncertainty in the error budget since we 
can use quite firm theoretical constraints on such evolution- 
ary parameter. As a fact, HB models available in the recent 
literature supply slightly discordant luminosity values, but 
very similar predictions concerning the mass. 

As well known, fundamental (RRab) and first-overtone 
(RRc) RR Lyrae populate the red and the blue side of the 
pulsation region, respectively. Therefore one may assume for 
F and FO pulsators the predicted stellar mass listed in Table 
1 for ZAHB models at logT e = 3.80 and 3.85, respectively. 
The total uncertainty introduced by this assumption is of the 
order of 3.5% if we simultaneously account for the predicted 
spread in mass between F and FO pulsators and for the 
^/^ccurrence of post-ZAHB evolutionary effects. 

By inserting these mass values into the previous re- 
log M^fM^ ns ' one eventually derives the following near-infrared 
Period- Luminosity-Metallicity (PLZk) relations: 

= 0.565 - 2.101 log P + 0.125 log Z- 0.734 log L (3) 

Mk° = -0.016 -2.265 log P + 0.096 log Z- 0.635 log L (4) 

with a rms scatter of <jk = 0.031 mag and 0.025 mag, re- 
spectively, including in quadrature an uncertainty of ~ 0.03 
mag, as given by the adopted spread in mass (3.5%) at fixed 
metal abundance. 

By following a similar procedure, one finds that also 



© 0000 RAS, MNRAS 000, 000-000 



4 G. Bono et al. 




Figure 1. Top panel: predicted My — Mk colors of F (filled 
circles) and FO pulsators (open circles) as a function of period 
(data in Table 2 and Table 3). Bottom panel: same as the top 
panel but with individual colors scaled to the same metallicity 
(logZ = —3), luminosity (log L/ Lq = 1.61), and mass (0.65Mq). 
The solid lines display the predicted relations (see equation (5) 
and (6) in the text). 



the T e dependence of the pulsator color can be replaced by 
a period dependence. In particular, since the V magnitudes 
are strongly correlated with L, whereas the K magnitudes 
do not, we expect that the V — K color is strongly depen- 
dent on luminosity, thus providing useful constraints to this 
unknown parameter. 

As a fact, using the whole set of pulsating models, we 
derive: 



{Mv-MkY = 4.014+1.986 log P-0.134 log Z-1.662 log L+l. 

(5) 

(M v -M K ) FO = 5.195+2.518 log P-0.159 log Z-2.158 log L+ 

(6) 

with a rms scatter of 0.027 mag and 0.023 mag, respectively 
(see Fig. 1). 

By assuming the predicted masses at logT e =3.80 and 
3.85 for F and FO pulsators, respectively, the previous rela- 
tions yield the Period- Luminosity-Color relations (PLCvk) 
in the V, K bands for both F and FO pulsators, as given by: 

(M V -M K ) F = 3.963+1.986 log P-0.162 log Z- 1.662 log L 

(7) 

(M v -M K ) FO = 5.180+2.518 log P-0.168 log Z-2.158 log L 

(8) 



with a total intrinsic dispersion of 0.037 mag and 0.031 mag, 
respectively, which accounts for the adopted spread in mass 
(3.5%) at fixed metallicity. 

According to eqs. (7) and (8), one finds that accurate 
V, K photometric measurements of F and FO pulsators with 
known period and metal content can supply independent in- 
formation on the intrinsic luminosity of the variable. More- 
over, since E(V - K)=0.89A V (Cardelli et al. 1989), if the 
visual extinction along the line of sight is known with an 
accuracy of ± 0.03 mag, then the luminosity might be eval- 
uated with a (formal) accuracy of ~ 0.03 dex. 



3 THE CASE OF RR LYR ITSELF 

As already discussed by Bono et al. (2002, Paper II), 
the predicted PLZk relation for F pulsators can be used 
to provide for the prototype variable RR Lyr a "pul- 
sational" parallax that appears in close agreement with 
the absolute value recently measured by Benedict et al. 
(2002, hereafter B02) on the basis of new astrometric 
data collected with FGS 3, the interferometer on board 
of HST. In that paper, we assumed an intrinsic luminos- 
ity of RR Lyr in the range logL/L0 = 1.65-1.8O to cover 
current uncertainties of HB models. On this basis, we de- 
rived Mk~— 0.541+0.062 mag and a "pulsation" parallax 
Kp U is=3. 858+0. 131 mas, in close agreement with the HST 
absolute value 7r a (, s =3.82±0.20 mas. 

According to the results presented in the previous 
section, one may revisit the case of RR Lyr, avoiding 
any assumption on the RR Lyr bolometric luminosity. 
Adopting logP=— 0.2466 (Kazarovets, Samus, & Durlevich 
2001), #=6.54 mag (Fernley, Skillen, & Burki 1993), 
[Fe/H]=-1.39±0.15 (Fernley et al. 1998a; Beers et al. 2000) 
and V=7.784 mag (Hardie 1955), equation (7) gives logi = 
1.642 ± 0.024 + 0.535A V . Thus, for each adopted extinction 
correction, one derives directly from the observed V — K 
color the luminosity value to be inserted into equation (3), 
without using any evolutionary predictions. Once Mk is 
determined from equation (3), both the intrinsic distance 
modulus (/tto = K — G.llAv — Mk) and the absolute visual 
magnitude (My = V — Av — no) can be easily determined. 

Data listed in Table 4 show that the pulsation par- 
allax presents a mild dependence on the adopted extinc- 
g5g^j5gy\ga5jpction. In particular, the top panel in Fig. 2 shows 
that when moving from Ay = 0.01 to 0.22 mag, the pul- 
sation parallax (filled circles) decreases from 3.89 ± 0.07 
to 3.78 ± 0.07 mas. However, current pulsational estimates 
within the errors of the HST absolute paral- 
lax (dashed lines), and present a better formal accuracy 
(/jjr/n ~ 3% versus ~ 5%). According to this evidence the 
pulsation parallaxes resemble direct astrometric measure- 
ments that are independent of interstellar absorption. 

Concerning the RR Lyr predicted absolute magnitude, 
the middle and the bottom panel in Fig. 2 show that 
when moving from Av=0.01 to 0.22 mag Mk changes from 
-0.513 ± 0.041 mag to -0.595 ± 0.041 mag, while M v from 
0.722 ± 0.041 to 0.453 ± 0.041 mag, in close agreement with 
HST results (open circles), but with a better formal accu- 
racy (±0.041 versus ±0.114 mag). In other words, despite 
the substantial improvement in the accuracy of the RR Lyr 
trigonometric parallax provided by HST, when compared 
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Figure 2. Top panel: predicted pulsation parallax for RR Lyr 
(filled circles) as a function of the visual extinction. Solid and 
dashed lines show the HST direct measurement and the relative 
uncertainty (B02), respectively. . Middle panel: comparison be- 
tween absolute if-magnitudes of RR Lyr based on HST trigono- 
metric parallax (open circles) and on the pulsation approach 
(filled circles), as a function of the visual extinction. For clarity 
sake, HST data are shifted by —0.005 along the x-axis. Bottom 
panel: the same, but for the absolute visual magnitude. The arrow 
marks the weighted mean of current extinction corrections 



with previous measurements (see B02) , data plotted in Fig. 
2 show that a sound empirical determination of its absolute 
magnitude is still hampered by the intrinsic uncertainty on 
the HST measurement, even if the interstellar extinction to 
RR Lyr would be firmly known. 

According to B02, the mean extinction correction to RR 
Lyr, based on photometric measurements of the astrometric 
reference stars located close to this variable, is < Av >= 
0.07 ± 0.03 mag. However, the same authors also supply an 
alternative value < Ay >= 0.11 ± 0.10 mag. Spectroscopic 
measurements of the diffuse interstellar band at A = 5780A 
(E(B-V) = -0.01 mag) and of the Na D-lines (E(B-V) = 
0.12 mag) carried out by Clementini et al. (1995) do not 
allow us to settle the problem of the reddening towards RR 
Lyr. On the other hand, empirical relations connecting the 
intrinsic color of RR Lyrae stars to the B-band amplitude, 
the period, and the metallicity (Piersimoni, Bono, & Ripepi 
2002) or to the Fourier amplitudes of the V-bsmd light curve 
(Kovacs & Walker 2001), together with data listed in Table 
5, supply a mean reddening of E(B~ V) = 0.029T0.005 mag. 
Finally, we note that the color excess given by Burstein & 
Heiles (1978) and by Blanco (1992) is E(B-V) = 0.03T0.02 
mag, while F98 adopts E(B — V) = 0.06 ± 0.03 mag. As a 
whole, the weighted mean of extinction corrections available 



Table 4. The case of RR Lyr itself. For each adopted extinction 
correction Ay-, column 2 gives the luminosity (in solar units) de- 
rived using the equation (7). This luminosity is used in equation 
(3) to evaluate Mr (3), and in turn to estimate together with 
K = 6.54 mag the intrinsic distance modulus, MO, and the "pul- 
sational" parallax tt pu i s . The last two columns give the absolute 
visual magnitude, based on V = 7.784 mag as well as on the 
pulsational and on the HST absolute parallax (7r a (, s =3.82±0.20 
mas), respectively. 



(mag) 


logL(7) 
±0.024 


M K (3) 
(mag) 
± 0.041 


Mo 
(mag) 
± 0.041 


(mas) 
± 0.072 


lvl v 
(mag) 

± 0.041 


M» ST 

V 

(mag) 
± 0.114 


0.01 


1.648 


-0.513 


7.052 


3.888 


0.722 


0.684 


0.04 


1.664 


-0.524 


7.060 


3.872 


0.684 


0.654 


0.07 


1.680 


-0.536 


7.069 


3.857 


0.645 


0.624 


0.10 


1.696 


-0.548 


7.077 


3.842 


0.607 


0.594 


0.13 


1.712 


-0.560 


7.086 


3.827 


0.568 


0.564 


0.16 


1.728 


-0.572 


7.094 


3.812 


0.530 


0.534 


0.19 


1.744 


-0.583 


7.103 


3.797 


0.491 


0.504 


0.22 


1.760 


-0.595 


7.111 


3.783 


0.453 


0.474 



Table 5. Empirical data for RR Lyr based on the photometry 
presented by Hardie (1955). The Fourier amplitudes Ai, A2 and 
A3 were estimated by fitting the F-band light curve by means of 
a cosine series with 7 components. 



p (d) 


0.5668 


(V) (mag) 


7.784 


(B-V) (mag) 


0.371 


A B (mag) 


1.154 


A v (mag) 


0.892 


Ai 


0.3014 


A 2 


0.1493 


A3 


0.0920 



in the literature gives < Ay >= 0.087T0.034 mag (see arrow 
in Fig. 2). 

By adopting this mean reddening value, the HST as- 
trometric measurement gives Mv = 0.607 ± 0.114 (intrinsic 
HST uncertainty) ±0.034 mag (A v contribution) = 0.607 ± 
0.119 mag. On the other hand, the pulsational approach 
gives logL = 1.689±0.030, and thus M v = 0.624±0.041 (in- 
trinsic uncertainty on the predicted Mk magnitude) ±0.044 
mag (Av contribution) = 0.624 ± 0.060 mag. 

The final uncertainty on the predicted RR Lyr param- 
eters appears still too large to firmly discriminate among 
current theoretical or empirical constraints on RR Lyrae in- 
trinsic luminosities and absolute visual magnitude. However, 
the use of the PLZk relation over a large number of well- 
studied RR Lyrae stars, see next section, that cover a wide 
metallicity range will provide relevant information on both 
the slope and, provided that no systematic errors are affect- 
ing the pulsating models, the zero-point of the My(RR)- 
[Fe/H] relation. 



4 FIELD RR LYRAE STARS 

Following the procedure discussed in the previous section 
for RR Lyr itself, the intrinsic luminosity of individual RR 
Lyrae stars with well-determined K and V magnitudes, as 
well as metallicity and reddening, can be directly inferred 
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Table 6. Observed parameters of RR Lyrae stars for which accurate visual and near-infrared light curves are available. The last column 
lists the sources of the data. The luminosity given in column (7) was derived using equations (7) [RRab] and (8) [RRc] and is in solar 
units. This luminosity is used in equation (3) [RRab] and (4) [RRc] to derive the Mx magnitude. On the basis of apparent and absolute 
K magnitudes we evaluate the intrinsic distance modulus, and in turn the absolute visual magnitude. 



ID logP [Fe/H] E(B-V) V K logL M K My Ref. a 

(mag) (mag) (mag) (mag) 

±0.04 ±0.05 ±0.07 
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a References: 1) Jones et al. 1988a; 2) Jones et al. 1988b; 3) Jones 
et al. 1992; 4) Liu & Janes 1990a; 5) Skillcn et al. 1993; 6) 
Cacciari ct al. 1992; 7) Fernley et al. 1990; 8) Storm et al. 
1994a; 9) Liu & Janes 1990b; 10) Storm ct al. 1994b. 
6 See text for more details. 



from the measured V ~ K color, without any a priori as- 
sumption related with the predicted luminosity of horizon- 
tal branch models. Once the individual luminosity of RRab 
variables is derived from equation (7) (equation (8) for RRc 
variables, the predicted PLZk relations [see equations (3) 
and (4)] provide the absolute near-infrared magnitude. This 
supplies the intrinsic distance modulus, and in turn the ab- 
solute visual magnitude. This means that we can investigate 
both the Mx-[Fe/H]-logP and the M v -[FeR] relations. 



We searched for RR Lyrae stars for which accurate V 
and K light curves, as well as reliable evaluations of both 
reddening and metallicity, are available in the literature. A 
total of 39 stars were found, with 30 field and 9 GGC mem- 
bers. Metal abundances have been taken from the homoge- 
neous compilation by Fernley et al. (1998a). The sample is 
listed in Table 6 together with key parameters such as the lu- 
minosity estimated according to equation (7) or (8), the Mk 
provided by equation (3) or (4), and the ensuing My magni- 
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Figure 3. Estimated individual luminosities of RR Lyrae stars 
in our sample versus metal abundance (see Table 6). Filled and 
open circles show field RRab and RRc variables, while filled tri- 
angles indicate cluster RRab variables. The asterisk and the open 
triangle show RR Lyr itself and the peculiar RR Lyrae V49 in 47 
Tuc, respectively. The error bars display current uncertainties in 
metal abundance and in luminosity. Dotted and solid lines refer 
to old and recent ZAHB models, respectively, while the dashed 
line depicts the exhaustion of the central He-burning phase. See 
text for more details. 



tudes. The errors associated with luminosity, near-infrared, 
and visual absolute magnitudes account for the intrinsic dis- 
persion of the predicted relations, as well as for an average 
uncertainty of ±0.20 dex on the metallicity and of ±0.03 
mag on Av- 

Figure 3 shows the derived "pulsational" luminosity of 
RR Lyrae stars in our sample as a function of the metal 
content. Filled and open circles display field RRab and RRc, 
respectively. RRab stars in GGCs (M4, M5, and M92) arc 
plotted as filled triangles, while RR Lyr itself with an aster- 
isk. The open triangle marks the position of the RR Lyrae 
V9 in 47 Tuc which was neglected, since it is peculiar (Storm 
et al. 1994b). Together with individual RR Lyrae stars, Fig. 
3 also shows the predicted log Lrr— [Fe/H] relation accord- 
ing to the old ZAHB models by Castellani, Chieffi, & Pulonc 
(1989, 1991, hereinafter CCP, dotted line) as well as the im- 
proved computations presented by Cassisi et al. (1998, 1999, 
hereinafter C99, solid line). The dashed line displays the lo- 
cus of central He exhaustion according to C99. 

As a whole, pulsational luminosities appear in reason- 
able agreement with the result of recent evolutionary com- 
putations based on improved input physics. In fact, with the 
exception of V9 in 47 Tuc, they show a dispersion in lumi- 
nosity of the order of 0.10 dex (dashed line) above the theo- 
retical ZAHB, in agreement with the predictions of synthetic 
horizontal branch simulations (see. e.g. Caputo et al. 1993; 
Demarque et al. 2000) . Thus it seems that recent evolution- 



Figure 4. Pulsational absolute visual magnitudes of RR Lyrae 
variables in our sample as a function of metal content. The sym- 
bols are the same as in Fig. 3. 



ary computations and pulsational constraints are supporting 
each other. 

When passing to magnitudes, one has to recall that stel- 
lar evolution predicts the luminosity and, therefore, the mag- 
nitude of "static" stellar structures, whereas the present pul- 
sational scenario refers to magnitude- weighted mean magni- 
tudes. As discussed by Bono, Caputo & Stellingwerf (1995), 
variables that present large luminosity amplitudes are ex- 
pected to have mean visual magnitudes up to ~ 0.15 mag 
fainter than the equivalent static values. Bearing in mind 
this caveat, Fig. 4 shows the predicted (magnitude- weighted) 
absolute visual magnitudes of the variables in our sample as 
a function of metallicity. The comparison with the ZAHB 
theoretical predictions by Cassisi et al. (1998, 1999) is now 
less meaningful and only discloses the compatibility of pulsa- 
tional and evolutionary scenarios, perhaps suggesting that 
the evolutionary models are somehow too bright with re- 
spect to pulsational magnitudes. However, a more detailed 
investigation concerning the comparison with current evo- 
lutionary predictions requires reliable amplitude corrections 
to the observed mean magnitudes to obtain the equivalent 
static values. 

It is worth noting that a least-squares solution through 
the data plotted in Fig. 4 suggests that the My- [Fe/H] rela- 
tion becomes steeper and steeper when moving from metal- 
poor to metal-rich variables (see Table 7). This result is in 
agreement with the change in the slope at [Fe/H]~ —1.5 
originally predicted by CCP and further supported by syn- 
thetic HB simulations (Caputo et al. 1993; Demarque et al. 
2000), by RR Lyrae stars in the field (MacNamara 1999), in 
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Table 7. The slope of the Mv(RR)-[Fc/H] relation as a function 
of the adopted metallicity range. 



[Fe/H] 


range 


slope 


-2.4/ 


-1.8 


»J 


-2.4/ 


-1.7 


0.080±0.100 


-2.4/ 


-1.6 


0.190±0.090 


-2.4/ 


-1.5 


0.192±0.079 


-2.4/ 


-1.4 


0.221±0.060 


-2.4/ 


-1.3 


0.282±0.053 


-2.4/ 


-1.0 


0.284±0.050 


-2.4/ 


-0.6 


0.283±0.039 


-2.4/ 


-0.3 


0.303±0.022 



GGCs (Caputo et al. 2000), and in to Centauri (Rey et al. 
2000). 



5 PLZ K VERSUS FOBE AND BW ABSOLUTE 
MAGNITUDES 

Figure 5 shows the comparison between the absolute visual 
magnitudes of field RR Lyrae stars estimated using the pre- 
dicted PLZ K relation (circles, see column (9) in Table 6) 
and the averaged absolute visual magnitudes of RR Lyrae 
stars in GGCs (triangles) derived by Caputo et al. (2000) ac- 
cording to the "First Overtone Blue Edge" (FOBE) method. 
The agreement between these two independent approaches 
is excellent, thus supporting the consistency of the pulsa- 
tional scenario. At the same time,this result indicates that 
field and cluster RR Lyrae stars should have experienced 
quite similar evolutionary histories. 

As a whole, linear interpolation through the data plot- 
ted in Fig. 5 gives the following analytical relations (solid 
lines) : 

M v = 0.718(±0.072) + 0.177(±0.069)[Fe/H] (9) 
for RR Lyrae stars with [Fe/H]< —1.6, and 

M v = 1.038(±0.077) + 0.359(±0.027)[Fe/H] (10) 

for RR Lyrae stars with [Fe/H]> —1.6. 

Let us now compare our results with current distance 
determinations based on the BW method. To perform a ro- 
bust comparison we adopted the absolute visual magnitudes 
for field RR Lyrae stars derived by Fernley (1994). The key 
feature of this sample is that the empirical data available 
in the literature were reduced by the previous author to the 
same NIR BW method. At the same time, for these objects 
the metal abundances were also estimated by adopting a ho- 
mogeneous metallicity scale (Fernley et al. 1998a). We ended 
up with a sample of 27 RR Lyrae stars. The nine GGC vari- 
ables were not included, since they appear to be affected 
by larger intrinsic errors and also because they have not 
been reduced to the same BW method yet. Figure 6 shows 
the comparison between these NIR BW absolute magnitudes 
(filled triangles) and the absolute visual magnitudes for the 
same objects, as based on the PLZk relation (circles). Data 
plotted in this figure display quite clearly that for RR Lyrae 
stars more metal-poor than [Fe/H] ~ —1.5, the discrepancy 




[Fe/H] 



Figure 5. Comparison between absolute visual magnitudes of RR 
Lyrae variables in our sample evaluated using the current pulsa- 
tion approach (circles) and the results of First Overtone Blue 
Edge (FOBE) method (triangles) to globular cluster variables 
(Caputo et al. 2000). The solid lines are the relations described 
in the text. 



is on average of the order of 0.2 mag, in the sense that BW 
magnitudes are fainter than those provided by current pul- 
sational approach. The discrepancy decreases toward higher 
metal contents and vanishes in the metal-rich tail. This in- 
dicates that the metallicity dependence given by the BW 
method is milder than predicted by pulsation models. 

Note that a new calibration of the BW method has been 
provided by Fernley et al. (1998b), but the difference over 
the entire metallicity range is negligible when compared with 
Fernley (1994). On the contrary, the preliminary analysis 
provided by Cacciari et al. (2000) of the assumptions cur- 
rently adopted in the BW method disclosed quite different 
results. In particular, by using available photometric data 
together with a large set of new atmosphere models, they 
found that RR Cet ([Fe/H] = -1.45) is « 0.12 mag brighter 
than previously estimated, whereas no significant correction 
was required for SW And ([Fe/H]=— 0.24). As a result, these 
new absolute magnitudes (open squares in Fig. 6), once con- 
firmed, appear fully consistent with the My — [Fe/H] rela- 
tion based on pulsation models [equation (10)]. This sug- 
gests that the slope of this relation, at least in the metal- 
licity range —1.5 <[Fe/H]< —0.4, is steeper than currently 
adopted (see also McNamara 1999). 

Following a referee's suggestion, one may investigate 
whether the discrepancy at [Fe/H] < — 1.5 between pul- 
sational absolute visual magnitudes and BW magnitudes 
might be settled by changing the mass value in equations 
(1) and (5). We found that to remove the difference of 
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-2 -1 

[Fe/H] 

Figure 6. Comparison between the absolute visual magnitudes 
estimated using the current pulsation approach and the indepen- 
dent evaluations derived by Fernley (1994) using the NIR BW 
method (triangles). The two open squares mark the position of 
RR Cet ([Fe/H]=-1.45) and SW And ([Fc/H]=-0.24) according 
to the new calibration of the BW method provided by Cacciari et 
al. (2000). The solid lines are the relations described in the text. 

~ 0.25 mag between BW and pulsational magnitudes at 
[Fe/H] = — 1.8, one should account for a decrease in mass of 
the order of ~ 23%, thus predicting a quite unrealistic low- 
mass, O.55M , for these metal-poor (Z ~ 0.0003) horizontal 
branch stellar structures. 

As for Jf-magnitudes, Fig. 7 shows that the absolute 
if-magnitudes estimated using the pulsational approach do 
obey a very tight relation with the metal content, once the 
dependence on period is removed. Owing to the small num- 
ber (3) of RRc variables, their periods have been fundamen- 
talized, (logJV=logJVo +0.127), thus yielding for the entire 
sample the following relation: 

M£ = -0.770(±0.044) -2.101 log P + 0.231 (±0.012) [Fe/H] 

(11) 

The small errors affecting the coefficients of this relation 
further 

strengthen the evidence that near-infrared observations 
could supply accurate distance determinations. 

To investigate in more detail the discrepancy between 
the current approach and the BW method, we compare the 
pulsational absolute K magnitudes with the absolute K 
magnitude obtained using the BW true distance modulus 
and the dereddened apparent K magnitudes listed in col- 
umn (6) of Table 6. Figure 8 shows the comparison between 
pulsational Mk values (circles) and results based on the 
BW method (triangles). Once again, the BW magnitudes at 
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Figure 7. Absolute X-magnitudcs of RR Lyrac in our sample 
estimated using the pulsation method as a function of metal con- 
tent. Data were projected onto a two dimensional plane. The solid 
line shows the projection onto the same plane of equation (11). 
The symbols are the same as in Fig. 3. 

[Fe/H]< — 1.5 are fainter than predicted by the pulsational 
approach, leading to a milder dependence on [Fe/H] with 
respect to the PLZk predictions (solid line). However, the 
absolute ^-magnitudes derived by Cacciari et al. (2000) for 
SW And and RR Cet (open squares) appear once again in 
good agreement with our pulsational approach. 



6 SUMMARY AND FINAL REMARKS 

Data plotted in Fig. 6 and Fig. 8 clearly show that the 
pulsational approach provides absolute magnitudes for RR 
Lyrae stars that are significantly brighter than suggested 
by the BW method. This difference is more evident for 
metal-poor ([Fe/H]< — 1.5) RR Lyrae variables. The re- 
cent revision provided by Cacciari et al. (2000) for RR Cet 
([Fe/H] = -1.45) and SW And ([Fe/H]=-0.24) agrees quite 
well with the pulsational results. This suggests that both 
the zero-point and the slope of the My- [Fe/H] and of the 
Mx-[Fe/H]-logP relation based on the old BW calibration 
need to be revisited. 

We have also found that both V and K absolute mag- 
nitudes of RR Lyr itself, estimated via the pulsational ap- 
proach, are in good agreement with the trigonometric par- 
allax recently measured by HST (see Fig. 2). This suggests 
that at least for the metallicity of RR Lyr ([Fe/H] = -1.39), 
the pulsational approach is consistent with direct distance 
determinations. 

As far as the metal- dependence is concerned, we finally 
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Figure 8. Comparison between the absolute K-magnitudcs based 
on the pulsation approach and the BW method. The latter were 
derived using the true distance moduli given by Fernley (1994) 
and the _ft"o- ma g n itudes listed in Table 6. The solid line refers to 
equation (11). 

note that the effect of metallicity on equation (11) can be 
revealed by a simple cross-check between predicted K mag- 
nitudes and periods. It has been already shown (see Paper 
I) that K measurements for RR Lyrae stars in M3, namely 
for variables at quite constant luminosity and metal con- 
tent, closely follow the pulsational constraints on the period 
dependence. On the contrary, data plotted in Fig. 9 show 
that for the current sample of cluster and field RR Lyrae 
stars the Mx-logP relation is significantly steeper than pre- 
dicted under the above quoted assumptions (dashed lines). 
This is the expected aftermath of the fact that metal-poor 
variables, which have longer periods (see data in Table 6), 
have also brighter bolometric luminosities (see Fig. 3) than 
metal-rich ones. 
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